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LAPLACIAN MATRICES AND SPANNING TREES OF TREE 

GRAPHS 

PHILIPPE BIANE AND GUILLAUME CHAPUY 


Abstract. If G is a strongly connected finite directed graph, the set TG of 
rooted directed spanning trees of G is naturally equipped with a structure of 
directed graph: there is a directed edge from any spanning tree to any other 
obtained by adding an outgoing edge at its root vertex and deleting the out¬ 
going edge of the endpoint. Any Schrodinger operator on G, for example the 
Laplacian, can be lifted canonically to TG. We show that the determinant 
of such a lifted Schrodinger operator admits a remarkable factorization into a 
product of determinants of the restrictions of Schrodinger operators on sub¬ 
graphs of G and we give a combinatorial description of the multiplicities using 
an exploration procedure of the graph. A similar factorization can be obtained 
from earlier ideas of C. Athaniasadis, but this leads to a different expression 
of the multiplicities, as signed sums on which the nonnegativity is not appear- 
ent. We also provide a description of the block structure associated with this 
factorization. 

As a simple illustration we reprove a formula of Bernard! enumerating span¬ 
ning forests of the hypercube, that is closely related to the graph of spanning 
trees of a bouquet. Several combinatorial questions are left open, such as 
giving a bijective interpretation of the results. 


1. Introduction 

Kirchoff’s matrix-tree theorem relates the number of spanning trees of a graph to 
the minors of its Laplacian matrix. It has a number of applications in enumerative 
combinatorics, including Cayley’s formula: 

(1.1) |rX„|=n"-\ 

counting rooted spanning trees of the complete graph with n vertices and Stan¬ 
ley’s formula: 

n 

(1.2) |r{0,iri = n(2*)(-t 

i=l 

for rooted spanning trees of the hypercube {0,1}") see [5]. In probability theory, a 
variant of Kirchoff’s theorem, known as the Markov chain tree theorem, expresses 
the invariant measure of a finite irreducible Markov chain in terms of spanning trees 
of its underlying graph (see P Chap 4] , or ( |2.3| ) below). An instructive proof of this 
result relies on lifting the Markov chain to a chain on the set of spanning trees of its 
underlying graph. In particular, this construction endows the set TG of spanning 
trees of any weighted directed graph G with a structure of weighted directed graph. 

Both authors acknowledge support from Ville de Paris, grant “Emergences 2013, Combinatoire 
a Paris”. G.C. acnowledges support from Agence Nationale de la Recherche, grant ANR 12-JS02- 
001-01 “Cartaplus”. 
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This construction is recalled in Section]^ (the reader can already have a look at the 
example of Figure]^. In the recent paper [3], the first author conjectured that the 
number of spanning trees of TG is given by a product of minors of the Laplacian 
matrix of the original graph G. In this paper, we prove this conjecture. More 
generally, given a Schrodinger operator on G, we will show (Theorem 3.5) that 
the determinant of a lifted Schrodinger operator on TG factorizes as a product of 
determinants of submatrices of the Schrodinger operator on G. In this factorization, 
only submatrices indexed by strongly connected subsets of vertices W C V{G) 
appear, and the multiplicity m{W) with which a given subset appears is described 
combinatorially via an algorithm of exploration of the graph G. 

The case of the adjacency matrix (another special case of Schrodinger operator) 
was already studied by C. Athanasiadis who related the eigenvalues in the graph 
and in the tree graph (see [5], or Section 3.1). As we shall see, this leads to a 
similar factorization of the characteristic polynomial as the one we obtain, and in 
fact the proof of [5] can easily be extended to any Schrodinger operator. However 
the methods of [2], whose proofs are based on a direct and elegant path-counting 
approach via inclusion-exclusion, lead to an expression of the multiplicities as signed 
sums which are not apparently positive. Our proof is of a different kind and proceeds 
by constructing sufficiently many invariant subspaces of the Laplacian matrix of 
TG. It is both algebraic and combinatorial in nature, but it leads to a positive 
description of the multiplicities. As a result our main theorem, or at least its 
main corollary, can be given a purely combinatorial formulation, which suggests 
the existence of a purely combinatorial proof. This is left as an open problem. 
Another combinatorial problem that we leave open concerns the definition of the 
multiplicities m{W): in the way we define them, these numbers depend both on a 
total ordering of the vertex set of the graph, and on the choice of a “base point” 
in each subset W, but it follows from the algebraic part of the proof that they 
actually do not depend on these choices. This property is mysterious to us and 
a direct combinatorial understanding of it would probably shed some light on the 
previous question. 

Finally, we note that there exists a factorization for the Laplacian matrix of the 
line graph associated to a directed graph (see [5]) that looks similar to what we 
obtain here for the tree graph. The case of the tree graph is actually more involved. 


The paper is structured as follows. In Section we state basic definitions 
and recall the construction of the tree graph. We also present the results of 
Athanasiadis [2] and rephrase them from the viewpoint of the characteristic polyno¬ 
mial. Then in Section we introduce the algorithm that defines the multiplicities 
m(W), which enables us to state our main result for the Schrodinger operators 
(Theorem |3.5[ ). We also state a corollary (Theorem |3.6| ) that deals with spanning 
trees of the tree graph TG, thus answering directly the question of [3]. In Sec¬ 
tion]^ we give the proof of the main result, that works, first, by constructing some 
invariant subspaces of the Schrodinger operator of TG, then by checking that we 
have constructed sufficiently enough of them using a degree argument. Finally in 
Section we illustrate our results by treating a few examples explicitly. 


Acknowledgements. When the first version of this paper was made public, we 
were not aware of the reference [2]. We thank Christos Athanasiadis for drawing 
our attention to it. G.C. also thanks Olivier Bernardi for an interesting discussion 
related to the reference [3]. 
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Figure 1. A directed graph G with 4 vertices (Left), and the 
graph TG (Right). Each vertex of TG is one of the 14 spanning 
trees of G. The weight of an edge in TG only depends on the 
root vertex of the two trees it links - only certain edge weights 
are indicated on this picture. The subset of vertices indicated next 
to each spanning tree is the i/i-value returned by the algorithm of 
Section 


3.2 


2. Directed graphs and tree graphs 

In this section we set notations and recall a few basic facts. 

2.1. Directed graphs. In this paper all directed graphs are finite and simple. Let 
G = (T, V) be a directed graph, with vertex set V and edge set E. For each edge 
we denote s(e) its source and t(e) its target. The graph G is strongly connected if 
for any pair of vertices {v,w) there exists an oriented path from v to w. 

If IF C F then the graph G induces a graph Gw = { W , Ew ) where Ew is the 
set of edges e with s{e),t{e) £ W. A subset W C V will be said to be strongly 
connected if the graph Gw is strongly connected. A cycle in G is a path which 
starts and ends at the same vertex. The cycle is simple if each vertex and each 
edge in the cycle is traversed exactly once. 

2.2. Laplacian matrix and Schrddinger operators. For a finite directed graph 
G, let Xe,e S T be a set of indeterminates. The edge-weighted Laplacian of the 
graph is the matrix (Qvw)v,wev given by Qyw = Xe if v ^ w s(e) = v,t(e) = w 
(this quantity is 0 if there is no such edge) and Qyy = — '^e-s{e)=v 

Let yy,v £ V he another set of variables and Y be the diagonal matrix with 
Yyy = yy. Tho associated Schrddinger operator with potential Y is the matrix 
L = Q -\-Y. Observe that, if one specializes the variables yy to a common value —z, 
then L = Q — zl and det(L) is the characteristic polynomial of Q evaluated on z. 

We will consider the space of functions on V with values in the field of rational 
fractions F^j = C{xe', e £ E,yy\v £V), and the space of measures on V (again with 
with values in F^j). These are vector spaces over the field F^j. The Schrodinger 
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operator L acts on functions on the right by 

L(j){v) = y^Ly^ip(w), 

W 

and on measures on the left by 

fJ.L{w) ='^fJ.{v)Lyyy. 

V 

The space of measures has a basis given by the 6v,v € V where Sy is the measure 
putting mass 1 on u and 0 elsewhere. 

2.3. A Markov chain. If the are positive real numbers, the matrix Q is the 

generator of a continuous time Markov chain on V, with semigroup of probability 
transitions given by This chain is irreducible if and only if the graph G is 

strongly connected. The function 1 is in the kernel of the action of Q on functions, 
and this kernel is one-dimensional if and only if the chain is irreducible. Dually, if 
the chain is irreducible then there is a positive measure in the kernel of the action 
of Q on measures (by the Perron-Frobenius theorem), which is unique up to a 
multiplicative constant. See for example [S] for more on these classical results. 

2.4. Spanning trees. Let G be a directed graph, an oriented spanning tree of G 
(or spanning tree of G for short) is a subgraph of G, containing all vertices, with no 
cycle, in which one vertex, called the root, has outdegree 0 and the other vertices 
have outdegree 1. If a is such a tree, with edge set i?a, we denote 

(2.1) ^ n 

More generally, if VF C F is a nonempty subset, an oriented forest of G, rooted in 
W, is a subgraph of G, containing all vertices, with no cycle and such that vertices 
in W have outdegree 0 while the other vertices have outdegree 1. Again for a 
forest f, with edge set Ef, we put 

(2.2) ^ n 

eGEf 

The matrix-tree theorem states that, if W C V and is the matrix obtained 
from Q by deleting rows and columns indexed by elements of W, then 

det((5'^) = TTf 

the sum being over oriented forests rooted in W. In particular, in the Markov chain 
interpretation, an explicit formula for an invariant measure is given by 

(2.3) fj.(v) = 

aeT„ 

where the sum is over spanning oriented trees rooted at v. This statement is 
known, in the context of probability theory, as the Markov Chain Tree theorem, 
see [HI Chap. 4]. 

It will be convenient in the following to use the notation Qvu = and 

Liy = denote the matrix extracted from the Laplacian or Schrodinger 

matrix of G by keeping only lines and columns indexed by elements of W. 
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2.5. The tree graph TG. Let G = {E,V) be a finite directed graph and a an 
oriented spanning tree of G with root r. For an edge e € V with s(e) = r, let b be 
the subgraph of G obtained by adding edge e to a then deleting the edge coming 
out of t{e) in a. See Figure]^ It is easy to check that b is an oriented spanning 
tree of G, with root t{e). 



Figure 2. An edge a — >■ b in the tree graph TG. It is associated 
with the edge weight Xe- 


The tree graph of G, denoted TG, is the directed graph whose vertices are 
the oriented spanning trees of G and whose edges are obtained by the previous 
construction, i.e. for each pair a, b as above we obtain an edge of TG with source 
a and target b. We will denote TV the set of vertices of TG, in other words, TV 
is the set of oriented spanning trees of G. Figure gives a full example of the 
construction. One can prove that the graph TG is strongly connected if G is, see 
for example [T]. Moreover the graph TG is simple and has no loop. There is a 
natural map p from TG to G which maps each vertex of TG, which is an oriented 
spanning tree of G, to its root, and maps each edge of TG to the edge e of G used 
for its construction. 

We assign weights to the edges and vertices of TG as follows: we give the weight 
Xe to any edge e' of TG such that p(e') = e and we give the weight to the tree 
a if its root is v. 

This leads to a weighted Laplacian and a Schrodinger operator for TG, which 
we denote respectively by Q and C. More precisely, Q is the matrix with rows and 
columns indexed by the oriented spanning trees of G such that 

Qac =0 if a T c and ac is not an edge of TG 

Qab =Xe if ab is an edge of TG and e is the edge of b going out the root of a. 

Qaa ~ ^ ^ Qah- 

b^a 

Similarly, y is the diagonal matrix indexed by TV with Taa = yroot(a) and 

c = Q + y. 

See [I] or [5] for more on the matrix Q in a context of probability theory. In [3] 
the first author proved that there exists a polynomial $(3 in the variables Xe such 
that, for any oriented spanning tree a of G, one has 

(2.4) det(Q‘^) = TTa^'G- 


In the same reference it was conjectured that $g is a product of symmetric minors 
of the matrix Q (i.e. a product of polynomials of the form det((5vF))- la this paper 
we prove this conjecture and provide an explicit formula for <i>G (Theorem 3.6). 
Actually we deduce this from a more general result which computes the determinant 
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of £ as a product of determinants of the matrices Lw (Theorem |3.5[ ). These results 
will be stated in Section and proved in Section The example of the tree graph 
of a cycle graph was investigated in [3] and we will explain in Section how it 
follows from our general resnlt. 


2.6. Structure of the tree graph. Before we state and prove the main theorem 
of this paper, we give here some elementary properties of the tree graph, which 
might be of independent interest. These properties will not be used in the rest of 
the paper. 

We start with the following simple observation: for any directed path tt in the 
graph G, starting at some vertex u, and any oriented spanning tree a rooted at n, 
there exists a unique path starting at a in TG which projects onto tt. Thus the 
graph TG is a covering graph of G. 

If a b is an edge of TG, then the union of the edges of a and b is a graph 
with a simple cycle G, containing the roots of a and b, and a forest, with edges 
disjoint from the edges of G, rooted on the vertices of G. The cycle G is the union 
of the path from the root of b to the root of a in the tree a with the edge from the 
root of a to the root of b in b. If we lift the cycle G in G to a path TG in TG, 
starting from a, we get a cycle in TG, which projects bijectively onto the simple 
cycle G. The cycle G, and thus TG is completely determined by the edge ab in 
TG, moreover for any edge in TG, the associated cycle is again TG. Conversely, if 
G is a simple cycle of G, and f a forest rooted at the vertices of G, then the trees 
obtained from G U f by deleting an edge of G form a simple cycle in TG which lies 
above G. We deduce: 

Proposition 2.1. The set of edges ofTG can partitionned into edge-disjoint simple 
cycles, which project onto simple cycles ofG. IfG is a simple cycle ofG, with vertex 
set W, then the number of simple cycles ofTG lying above G is egual to the number 
of forests rooted in W. 

In particular, to any outgoing edge of a in TG one can associate the incoming 
edge of the cycle to which it belongs, and this gives a bijection between incoming 
and outgoing edges of a. An immediate corollary is 

Corollary 2.2. The graph TG is Eulerian: The number of outgoing or incoming 
edges of a vertex a are both egual to the number of outgoing edges of the root of a 
in G. 


The previous discussion also implies that the measure on vertices of TG that 
gives mass TTa to each tree a is an invariant measure on TG, i.e. one has nR — 0. 
By using the projection map p, it follows that the measure fi on V given by (2.3) 
satisfies p,Q = 0, which gives a simple proof of the Markov Chain tree theorem, see 

m- 


3. A FORMULA FOR THE DETERMINANT OF THE SCHRODINGER OPERATOR 

We use the same notation as in the previous sections, in particular V is the 
vertex set of the directed graph G, the weighted Laplacian of G is Q, its Schrddinger 
operator is L, the graph of spanning trees is denoted TG and the weighted Laplacian 
and Schrddinger operators of TG, as in Section |T5| are denoted by Q and £. We 
assume that G is strongly connected. 
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3.1. Eigenvalues of the adjacency matrix, according to Athanasiadis [2]. If 

the weights yy are set to yy = —Qyy = J2e-s(e)=v then the Schrodinger operator 
becomes the adjacency matrix of the graph G. We denote it by M. It is easy to 
see that in this case the lifted Schrodinger operator Ai is the adjacency matrix of 
the graph TG. In [2] C. Athanasiadis proves the following result about eigenvalues 
of the matrix AA. 


Proposition 3.1 ([2]). The eigenvalues of the adjacency matrix AA are eigenvalues 
of the matrices Mx',X C V. For such an eigenvalue 7 , if mxij) denotes its 
multiplicity in Mx, then its multiplicity in AA is 

TUX (7) det(rx - I) 

XCV 

where Fx is the matrix Mx with all variables Xe equal to -1. 


The previous theorem implies the following equation 

det(z/-M)= Yl det{zl - MxY^^^ 

XCV 

where 1{X) = det(rx — !)■ Observe however that the multiplicities 1{X) can be 
negative in this equation. In order to get nonnegative multiplicities, we will use 
the following fact which is easy to check: for any A C P, if we let X = UiWi be 
its decomposition into strongly connected components, then the graph induces an 
order relation between the Wi from which one deduces the factorization 

det( 2 ;/ — Mx) = det(z/ — Mwi)- 

i 

It follows that 

(3.1) det{zl - AA) =Y\dei[zI - 

w 


where the product is over strongly connected subsets W <ZV and 
(3.2) n{W)= Y. 




where A Dsc IF means tha t IF is a strongly connected component of A. As 
we will see later (Lemma |4.l[ ), the polynomials det( 2 ;/ — Mvp)"*-^^ for IF strongly 
connected are distinct prime polynomials, therefore the formula |3.1| uniquely defines 
the multiplicities n{X) which therefore are nonnegative integers. This property 


however is not apparent from the formula 3.2 

In this paper, we will generalize this result to the case of Schrodinger operators 
and give another expression for the multiplicities, as the cardinality of a set of com¬ 
binatorial objects (hence the nonnegativity will be apparent). We will also explicitly 
exhibit a block decomposition of the matrix C that underlies the factorization of 
the characteristic polynomial. 

Although Athanasiadis’s results were stated for adjacency matrices, his proof 
actually extends easily to the more general case of Schrodinger operators which we 
consider here (with the same multiplicities). However the link between the two 
approaches is yet to be understood. 
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3.2. The exploration algorithm. Our formula for the determinant of £ (given in 
Theorem 3.51 involves certain combinatorial quantities dehned through an algorith¬ 
mic exploration of the graph. The exploration algorithm associates to any spanning 
tree a of G two subsets of vertices of G, denoted by (p^a) and i/’(a). Roughly speak¬ 
ing, the algorithm performs a breadth first search on the graph G, but only the 
vertices that are discovered along edges belonging to the tree a are considered as 
explored. Vertices discovered along edges not in a are immediately “erased”. This 
may prevent the algorithm from exploring the whole vertex set and, at the end, 
we call (j){a) the set of explored vertices. The set tjj{a) is the strongly connected 
component of the root vertex in (j){a). 

We now describe more precisely the algorithm. Because it is based on breadth 
first search, our algorithm depends on an ordering of the vertices of V. This ordering 
can be arbitrary but it is important to fix it once and for all: 


From now on we fix a total ordering of the vertex set V of G. 

In particular on examples and special cases considered in the paper, if the vertex 
set is an integer interval, we will equip it with the natural ordering on integers 
without further notice (this is the case for example on Figure [^. 


Exploration algorithm. 

Input: A spanning tree a of the directed graph G = {V, E), rooted at v. 

Output: A subset of vertices ^(a) C V; 

A subset of vertices tpia) C </>(a), such that G|^(a) is strongly connected. 

Running variables: - a set A of vertices of G; 

- an ordered list L of edges of G (first in, hrst out); 

- a set F of edges of G. 

Initialization: Set A := {u}, F := {e G E\s{e) u}, and let L be the list of edges 
of G with target w, ordered by increasing source. 

Iteration: While L is not empty, pick the hrst edge e in L and let w be its source: 

If e belongs to the tree a: 
add w to A; 

delete all edges with source w from L; 

append at the end of L all the edges in F with target w, by increasing source, 
else 

delete from L and F all the edges with source or target w in E. 

(in this case we say that the vertex w has been erased) 

Termination: We let (j){a) := A be the terminal value of the evolving set A. The 
directed graph G induces a directed graph on ^(a), and we let ^(a) be the strongly 
connected component of v in this graph. 


Observe that if a vertex w is picked up by the algorithm at some iteration, it 
will not appear again, this implies that the algorithm always stops after a finite 
number of steps. We refer the reader to Figure for an example of application 
of the algorithm. The reader can also look back at Figure on which, for each 
spanning tree a, the value of the set 'tp{a) is indicated. 

With the exploration algorithm, we can now define the multiplicities that are 
necessary to state our main theorem. 
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^ = { 3 } 




^(a) = {3,4} 


Figure 3. Left: in plain edges, a spanning tree a of the graph 
of Figure We initialize the set A to {3} since 3 is the root 
of a. Center: at the first two iterations of the main loop of the 
algorithm, we consider the edges (2,3) and (4,3), that belong to 
the tree a: the vertices 2 and 4 are thus added to the set A. Right: 
at the next iteration, we consider the edge (1,2) that does not 
belong to a. The vertex 1 is thus erased. The set A will not evolve 
until the termination step, and we thus get (j>{a) = {2,3,4}. The 
strongly connected component of 3 inside {2,3,4} in the original 
graph is {3,4}, which gives the value of tp{a.). 


Definition 3.2. Let W be a strongly connected subset of V, and w G W. The 
multiplicity ofW at w is the number m{W, w) of oriented spanning trees a rooted 
at w such that ip{a) = W. 

For any v G V, there exists a unique tree ay^v rooted at v such that ip{av,v) = V. 
This tree is obtained by performing a breadth first search on G starting from v and 
keeping the edges of first discovery of each vertex. We thus have: 

Lemma 3.3. For any v gV one has m{V,v) = 1. 

More generally, we will prove in Section |4.5| the following fact 

Definition-Lemma 3.4. For any strongly connected subset W G V, the multiplic¬ 
ity miW,w) depends neither on w G W nor on the ordering of the elements ofV. 
We will call m{W) this common value. 

Proof. See Section [431 ^ 

3.3. Main result. Our main result is the following theorem. 

Theorem 3.5. Let G be a strongly connected directed graph. Then the determinant 
of the lifted Schrodinger operator on TG is given by: 


(3.3) det{C) = Yl detiLwr^^^ 

wcv 

Ws.c. 

where the product is over all strongly connected subsets W QV. 


From the previous result we will deduce the following formula for $< 3 . Recall 
that we defined TTa in (2.1 ) as the product over the weights of the edges of a tree a 
and similarly TTf (2.2) for a forest. Analogously one defines the weight of a spanning 
tree of TG as the product of the weights of its edges. We define the polynomials 
Pgj Ftg and 'I'ty as the sums of these weights over, respectively, spanning trees of 
G, of TG, and of forests of G rooted in W. The Markov chain tree theorem implies 
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that the generating function of the spanning trees of a graph is the coefficient of the 
term of degree 1 in the characteristic polynomial of the Laplacian matrix. Using 
this fact and Theorem |3. 5 1 we obtain the following result. 

Theorem 3.6 (Spanning trees of the tree graph). The generating polynomial F-fc 
of spanning trees of the tree graph is given hy 


(3.4) 

Ftg = ^gFg, 

where 


(3.5) 

= n (^nvu) 


wcw 

Ws.c. 


where the product is over all proper strongly connected subsets W C V. 


Note that from (2.4) and the matrix-tree theorem, Theorem 3.6 also gives a 


formula for spanning trees of TG rooted at a particular spanning tree a. 


Note also that summing over all trees a in (2.4) and using the matrix-tree theo¬ 


rem, we see that the constant in (2.4) is indeed the same as the one in (3.4). 


Remark 3.7. Both sides of Equation (3.5) have a natural combinatorial meaning; 


the left hand side is a generating function for spanning trees of TG, while the right 
hand side is the generating function of some tuples of forests on G. It would be 
interesting to have a direct combinatorial proof of this identity. 


As an example, on the graph G of Figure[2 there are 7 strongly connected proper 
subsets of vertices and we have: to({1}) = ?7 i({2 }) = to({3}) = 0, to({ 4}) = 3, 
m({3,4}) = 2, m({l,2,3}) = 0, to({ 1,3,4}) = 1, and miV) = 1. It follows that 
the characteristic polynomial of the Schrodinger operator of the graph TG in this 
case is given by 

det(£) = det(L 4 )^ det(L 3 , 4 )^ det(Li_ 3 ^ 4 ) det(L). 

This identity can of course also be checked by a direct computation. 


4. Proof of the main results 


In this section we prove the main results. We assume as above that G is strongly 
connected and we use the same notation as in previous sections. 


4.1. Polynomials. In order to prove Theorem |3.5| we will show that each factor 
in (3.3) appears with, at least, the wanted multiplicity and conclude by a degree 
argument. We start by showing that these factors are irreducible. 


Lemma 4.1. IfWcV is a proper strongly connected subset then the polynomial 
det(Lvu) *5 irreducible as a polynomial in the variables {xe)eGE', {yv)vew- 

Proof. First we note that det(Lw) is a homogeneous polynomial, and it has degree 
at most one in each of the variables Xe,e G E,s{e) G W,(yv)y^w- Moreover, by 
Kirchhoff’s theorem, its term of total degree 0 in the y variables is the generating 
function of forests rooted in U \ W, which is nonzero since W is strongly connected 
and proper. In particular, the polynomial is not divisible by any of the yy. By 
expanding the determinant det(Lvu) along the row indexed by some w G W, we see 
that for each w, in each monomial of det(Lvi/) there is at most one factor Xg with 
s(e) = w. It follows that, for each w, as a polynomial in the variables (xg) s(e) = w), 
the polynomial det(Lvu) has degree 1. 
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Now assume that det(LvF) = AB is a nontrivial factorization into homogeneous 
polynomials then, from the previous point, for each w the polynomial AB is a 
factorization of a degree one polynomial in {xe',s{e) = w). It follows that there 
must exist a partition of W = X' where A is a polynomial in the i/„ and in the 
variables Xe with s(e) S X, while i? is a polynomial in the and in the variables Xe 
with s(e) S X'] note that this partition is non trivial since det(Lw) is not divisible 
by any ?/„. Moreover every monomial of det(LvF) can be written in a unique way as 
a product of a monomial appearing in A and a monomial appearing in B. Putting 
all variables x^ with s(e) & X to zero we see that 

det(Lw)|^^^o = det(Lx') Wyv= A{y,0)B. 

vGX 

The same can be done for X' and we obtain that det(Lx) Aei{Lx') = h{y) det(Lw) 
where h{y) = A{y, 0)B{y, 0) n„eVF ^ is a Laurent polynomial. By looking at the 
top coefficient in the y^ on both sides it follows that h = 1, hence 

det{Lw) = det(Ljc) det(Lx')- 

Since the graph Gw is strongly connected there exists a spanning tree a of IT rooted 
in some vertex x G X; in the corresponding monomial term of det(Lw) there is a 
factor Xe with s(e) = x' for each x' G X', since each vertex of X' has an outgoing 
edge in the tree a. The corresponding monomial therefore appears in det(Lx')) ^cid 
we note that each variable Xe appearing in this monomial is such that t(e) € W. 
The argument can be repeated for X and we deduce that there exists a monomial 
term in det(Lw) = det(Lx) det(Lx') which is a product of variables Xe which are 
all such that t{e) G IT; this monomial does not correspond to a forest by a simple 
counting argument, hence a contradiction. □ 

4.2. The case of the full minor. The space of functions on TG which depend 
only on the root of the tree (i.e. functions F such that F{sl) = F{h) if p(a) = p(b)) 
is invariant by the action of C on functions, moreover the restriction of C to this 
subset is clearly equivalent to the action of L on the functions on V by the obvious 
map. Dually the matrix C leaves invariant the space of measures /r such that 
y{p~^{v)) = 0 for all v G V. The action of C on the quotient of meas{TG) by 
this subspace is isomorphic to the action of L on meas(G). From either of these 
remarks, we deduce 

Lemma 4.2. The polynomial det{L) divides det{C). 

4.3. Boundary and erased vertices. We make some remarks on the algorithm 
of Section [3.2[ Once we have applied the algorithm to a given tree b, with output 
IT = 'i/'(t>), we can distinguish several subsets of vertices: 

(1) the set Z = ^(b), which is the set of vertices of a subtree of b; 

(2) the set IT = V'(t>), which is the set of vertices of a subtree of the previous 
one; 

(3) the set T = V\Z\ 

(4) the set of erased points which are the vertices which have been erased when 
applying the algorithm. 

(5) the set of boundary points, which are the vertices in Y having an outgoing 
edge with target in Z. 

Lemma 4.3. The sets of boundary points and of erased points coincide. 
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Proof. In an iteration of the algorithm, any vertex which has been added to the set 
A has all its outgoing edges suppressed, therefore it cannot be erased in a subsequent 
iteration. It follows that, if a vertex has been erased during the algorithm, then it 
does not belong to Z and it is the source of some edge with target in Z therefore it is 
a boundary point. Conversely if u is a boundary point let z € Z be the first vertex, 
among the targets of an outgoing edge of v, which is scanned by the algorithm, 
then the edge from v to z does not belong to the tree b (if it did, v would be in Z), 
therefore v is erased when one applies the algorithm at z. □ 


4.4. Constructing the invariant subspaces. Let W C V he & strongly con¬ 
nected proper subset. In this section and the next we will construct m(W) comple¬ 
mentary vector spaces that are invariant by £ and on which £ acts as the matrix 
Lw- This will be the main step towards proving (3.3). This construction goes in 
two steps: we first build a space of measures that is not invariant (this section. 


4.4) and we then construct a quotient of this space by imposing suitable “boundary 
conditions” that make the quotient space invariant (Section |4.5[ ). 

For every pair (a, f) formed of a spanning oriented tree a of IF and an oriented 
forest f rooted in W, let us call a x f the oriented spanning tree of V, rooted in 
the root of a, obtained by taking the union of the edges of a and f. Let us denote 
by TW the set of oriented spanning trees of W and TW the set of oriented forests 
rooted in IF. We thus have an injection TIF x TW TV and correspondingly 
a linear map from meas{TW) 0 meas{TW) —>■ meas{TV). Fix some forest f as 
above and consider the matrix £(vu) obtained from £ by keeping only the rows and 
columns corresponding to oriented spanning trees of V of the form a x f where a is 
some spanning oriented tree of IF. It is easy to see that this matrix, considered as 
a matrix indexed by elements of TIF, does not depend on the forest f, but only on 
IF. It differs from the matrix £ constructed from the graph Gw by some diagonal 
terms corresponding to the fact that there exists edges in E with source in IF and 
target in F \ IF. The matrix £(w) acts on functions on TIF and on measures on 
TIF, and it is easy to see that for its action on measures, the space of measures 
on TIF such that = 0 for every vertex w G IF is an invariant subspace 

of measures. The action of £(vu) on the quotient of meas{TW) by this subspace is 
isomorphic to the action of Lw on meas(IF). 


4.5. Boundary conditions and proof of Theorem 3.5 

sures meas{TW) Z> meas{TW) 


The subspace of mea- 


C meas{TV) is not invariant by the action of £ 
on measures but we will see that by modifying it and imposing suitable ’’boundary 
conditions” we will obtain an invariant subspace. For this let us consider a vertex 
w e IF and a tree b, rooted at w, such that '(/’(b) = IF. The tree b is of the form 
a X f considered above, moreover the tree a depends only on IF and w, since it co¬ 
incides with the breadth-first search exploration tree on IF (similarly as in Lemma 
To emphasize this fact we use the notation a = S-w.-u 


3.3). 


The set of trees b 

rooted at w and such that '(/’(b) = IF is equal to aw,w x Fw,w where Fw,w is some 
set of forests rooted in IF, with |Tiv,u;| = m(W,w). As indicated by the notation, 
the set Fw,w may depend on both IF and w. 

Let us fix f S Fw,w and consider the set €b of vertices erased when running the 
algorithm on the tree b = aw,v x f. A vertex v is erased when it is the source of 
some edge e{v) considered in the algorithm, which is not in b and which is scanned 
before the edge of b going out of v. For a subset S’ C £b let % be the graph 
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obtained by replacing in f, for each erased vertex v G the edge going out of v by 
the edge e{v). 

Lemma 4.4. For each 5^ C l£b the graph f;j is a forest rooted in W 

Proof. It suffices to observe that each vertex oi V \ W has outdegree 1 and that, 
by construction, from any such vertex there is directed path going to W. □ 


For f S let Vf be the measure on PW defined by 

(4.1) z/f = ^ 

5C(£b 

with b = an/,t) x f. 

Lemma 4.5. The measures Vf for f in Fw,v o-re linearly independent. 


Proof. First, construct a gradation on the set of spanning trees of G rooted at w 
as follows. If a is such a tree, let Vq = w,Vi,.. .Vk be the list of elements of ^(a) 
(the set of non-erased vertices) in the order they are discovered by the algorithm 
running on a. We let Iq, .. . ,lk be the number of incoming edges of uq, ..., Ufc from 
the set (^(a). This construction associates to any tree a rooted at in a finite sequence 
Iq, ... ,lk of integers. We equip the set of all sequences with the lexicographic order, 
which induces a gradation on the set of trees rooted at w. 

Now, if f G and S' ^ 0, then the tree a.w,w x % is strictly higher in 

the gradation than b.w,w x f. Indeed the origin of the hrst edge e{v) for u G S 
that is considered by the algorithm belongs to the set (j){au],w x %) but not to 
4>{aw,w X f), which shows that at the first index where the degree sequences differ, 
the one corresponding to sl.ui,w x % takes a larger value - hence it is larger for the 
lexicographic order. 

This shows that the transformation (4.1) expressing the measures {pg, g G F\y,w} 
in the basis {(5f,f G PW'\ is given by a matrix of full rank: indeed, provided we 
order rows and columns by any total ordering of PW that extends the gradation 
defined by f < g if {siw,w x f) < {a.\Y,w x g), we obtain a strict upper staircase 
matrix. □ 


It follows from the last lemma that the collection of measures 

i5t ® Pf 

where t runs over all rooted spanning trees of W and f over all elements of Fw,w 
is a linearly independent family of measures on PG. 

Now fix as above a forest f G Fw.w and let b = aw.w x f. Recall that = W, 
and call B C V \ W the set of boundary points relative to the tree b, as dehned in 
Section |4.3[ Let FI be the subgraph of TG where we have erased all edges having 
for source a tree rooted in a vertex of B. Let K be the subset of vertices of H 
which can be reached by a path in H starting from a tree of the form t x f;j, for 
some spanning tree t of IF and some 5^ C £b- Let J <Z K he the subset of trees 
whose root is not an element of B. Note that B, H, K and J all depend on the 
choice of f (or b) even though we do not indicate it in the notation. 

Lemma 4.6. Let £f be the space of measures spanned by dt 0 v^CP for all spanning 
trees t of W and all n> 0, then every measure in this £f is supported on the set J. 
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Proof. It is enough to prove that for all t and n > 0, the measure St 0 Uf£" has 
support in J, since this property is preserved under taking linear combinations. Let 
us compute St 0 Uf£"(c) for a tree c rooted a some boundary point v £ B. Recall 
that boundary points and erased points coincide by Lemma |4.3[ One has 

(4.2) <5tOuf/:”(c) = 

see IT 


where the sum is over all paths tt of length n in BG starting at t x fg and 

ending in c and £^ is the product of £e over all edges e traversed by tt. Let tt be 
such a path and r its projection on G, then the quantity is equal to Lr- Assume 
that V ^ S then the path r can be lifted to a path tt' starting at t x The 

only difference between f:g and /gu{D} is the edge coming out of v. Since the path 
T ends in v, the edge starting from v is deleted in the end tree of tt', therefore the 
end point of tt' is again c. It follows that the contributions of £^ and to the 
sum cancel. If u € 1? we consider the path tt' started at t x again the two 

contributions cancel. It follows that any contribution to the right hand side of (4.2) 
comes with another which cancels it, therefore the quantity (5t 0 Uf£"(c) vanishes 
for all n and all trees c rooted in some boundary point. 

Let now c be a tree which does not belong to the set K. We prove that 


(4.3) 


St 0 Uf£"(c) = 0 


by induction on n. Clearly this is true if n = 0 and 

dt <8 Uf£”+^(c) = (g) Uf/:"-(d)£dc. 

d 


Since c ^ K, if (£)dc 7 ^ 0 then either 

i) d is rooted in a boundary point, or 

ii) d ^ K. 

In the first case St 0 Uf£"’(d) = 0 by the first part of the proof. 

In case ii) St 0 Uf£”(d) = 0 follows from the induction hypothesis. 
Equation (4.3) follows. 


□ 


Now we let £ be the span of the spaces £f for all forests f S Fw,w Equivalently 
£ is the space of measures spanned by [St0vt\C'^ for all t, n and f. By construction 
the space £ is invariant by the action of £ on measures. 

Lemma 4.7. The subspace P of £ which consists of measures supported by trees 
with root not in W is an invariant subspace. 


Proof. It is enough to prove that for each f € Fw.w the subspace of £f which consists 
of measures supported by trees with root not in W is an invariant subspace. This is 
clear from the last lemma, since in the graph FI we have suppressed edges coming 
out from vertices of B (boundary vertices), hence it is not possible for a path to 
come back in W after having left it. □ 


Lemma 4.8. The action of C on the quotient space TfF carries m{W,w) copies 

of F(w) ■ 


Proof. Indeed for each forest f in F]y,w and any b spanning rooted tree of W, the 
measure 0 Uf satisfies [i5b 0 Uf]£ = [(Ib£(iv)] 0 + X where x £ Moreover 

the space span{vf) has dimension to(IE, re) by lemma 4.6 The lemma follows. □ 
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We can now finish the proof of the main results. 


Proof of Definition-Lemma \3.4\ and Theorem \3.5[ From Lemma [4~8| it follows that 
det(£) is divisible by 

for any strongly connected W and w G W. In particular we can take m{W,w) to 
be maximal among all w in W. This implies, since the different det{Lw) are prime 
polynomials (see Lemmas 4.11 that 

(4.4) 

is divisible by 

(4.5) 


det(£) 

det(L) X det(LvF)““”'^"'^’”^’^’“^^ 


w s.c. 


Now, the degree of (4.4) is \TV\ while that of (4.5) is 1^1 inaxu,gi^(m(IF, w)), 
therefore 

\TV\ > |IF| max(TO(IF, w)). 

w^W 

W 

By definition of m(IF, w) we have: 

I'^^i = E E 

w wew 

It follows that we have equality for all w € IF: 

m{W,w) = max(m(IF, li;)). 

w^W 

This proves that m(IF, w) does not depend on ru S IF and justifies the notation 
m(IF). This also proves that m(IF) is the multiplicity of the prime factor det(Lvy) 
in det(£). This quantity does not depend on the order chosen on F, thus justifying 
Definition-Lemma 13.41 


We have thus proved that the two sides of (3.3) are scalar multiples of each 
other. The proportionality constant is easily seen to be 1 by looking at the top 
degree coefficient in the variables y. □ 


5. The case of multiple edges 

Although Theorem |3.5| only covers the case of simple directed graphs, it is easy 
to use it to address the case of multiple edges. Indeed there is a well-known trick 
which produces a directed graph with no multiple edges, starting from an arbitrary 
directed graph, which consists in adding a vertex in the middle of each edge of 
the original graph. These new vertices have one incoming and one outgoing edge, 
obtained by splitting the original edge. This produces a new graph G = (F, E) 
with |F| = |F| + |if| and |if| = 2|i?|. Given a vertex v G V there is a natural 
bijection between spanning trees of G and of G rooted at v. For a vertex v of the 
new graph sitting on an edge e with s(e) = v oi G, there is a natural bijection with 
the spanning trees rooted at v. Thus the graph TG is obtained from TG by adding 
vertices in the middle of the edges. It is now an easy task to transfer results on G 
to results on G. We leave the details to the interested reader (the examples of the 
next section may serve as a guideline for this). 

Note that we do not need to take care of loops, that are irrelevant to the study 
of spanning trees. 
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6. Examples and applications 
In this section we illustrate our result on a few simple examples. 

6.1. The cycle graph. This example was treated in [3], let us see how to recover 
it via our main result. Let G = (V,E) be the cycle graph of size n, with vertex 
set V = [l..n] and a directed edge from itojiij = i±l mod n. Thus G has 
n vertices and 2n directed edges. The graph G has spanning trees: a spanning 
tree a is characterized by its root vertex r G [l..n] and by the unique i € [l..n] such 
that {i, i + 1} mod n are the two vertices of degree 1 in the tree. 

We note that for any subset of vertices W C V of cardinality n — 1, one has 
m{W) = 1. To see this, recall that m{W) = m{W,w) for any w GW and choose 
for w a neighbour of the unique vertex u not in W: then it is clear that the only 
spanning tree a such that '!/i(a) = IT is the one rooted at w in which u and w have 


degree 1. It is then easy to see, either directly or by considering the degree of (3.3), 
that these are the only proper subsets W C V such that m{W) ^ 0. 

Applying Theorem |3.6[ we obtain that is the product of all symmetric minors 
of Q of size n — 1, which was Theorem 2 in [4]. 

6.2. The complete graph (spanning trees of the graph of all Cayley trees). 

If G = Kn is the complete graph on n > 1 vertices, then TG is the set of all rooted 
Cayley trees of size n, thus TG has vertices by Cayley’s formula (1.1). If a is 


a Cayley tree rooted at r € [l..n], applying the exploration algorithm to a has the 
following effect: at the first step, all neighbours of r in a are explored and added 
to A, and all other vertices of T\ {r} are erased. It follows that for any W C [l..n] 
and w € W, the multiplicity m{W, w) is equal to the number of Cayley trees rooted 
at w in which the root has 1-neighbourhood W \ {w}. Those trees are in bijection 
with spanning forests of [l..n] \ {w} rooted at 1T\ {w}. We obtain, using a classical 
formula for the number of labeled forests of size n — 1 rooted at fc — 1 fixed roots: 

m{W) = m{w, W) = (k— l)(n — 1)"“^“^, where k = |1T|. 

This formula for the multiplicity appeared as a conjecture by the second au¬ 
thor in |3]. It is however easily seen to be equivalent to an earlier result of 
Athanasiadis [21 Corollary 3 . 2 ], which also refers to an earlier conjecture of Propp 
(we were not aware of the reference [2] at the time |1] was written). By applying 
Theorem 3.6 we obtain that the number of spanning trees of the graph TKn is 
equal to: 

n—1 


n {{n-ky 






k=l 


It would be interesting to give a direct combinatorial proof of this formula. 


6.3. Bouquets, and the hypercube. Fix k>l and integers ni,n 2 , ■ ■■ ,nk > 1. 
Consider the bouquet graph B with vertex set 

V = {0,1,..., A:} l±) , 1 < i < fc, 1 < j < nJ 

and a directed edge between each and 0, between each vertex i and each vj for 
1 < j < Hi and between 0 and each vertex i in [l..fc]. See the following picture: 
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For 1 < 7 < fc, 1 < j < rii, we assign the weight xl to the edge entering the 
vertex vl, the weight Si to the edge going from 0 to the vertex i, and we assign 
the weight 1 to all other edges. A spanning tree of B rooted at 0 is naturally 
parametrised by the index in of the edge outgoing from each vertex i in 

[l-./c]. We let am be the spanning tree rooted at 0 naturally parametrized by 
m € X [l..n 2 ] x • • • x [l..nk]- For each m, the tree has k outgoing edges 

in TB, to trees that we note 6^, for i £ where 6^ rooted at the vertex i 

and m' is the projection of m to x [l..n 2 ] x • • • x x • • • x [l..nfc] (i-th 

set in the product omitted). Each tree 6)^, has an outgoing path of length 2 going 
to each tree a^ such that m projects to m'. For example if A: = 1, then T~B is the 
following “star graph”: 



For k > 1, TB can be interpreted as a “partial product” of such star graphs of 
parameters ni,n 2 , ■ ■ ■ ,nk, more precisely it is the subgraph of the product of these 
graphs induced by the subset of vertices that are such that at most one of their 
coordinates is a vertex which is not “of type a”. In particular, if ni = n 2 = • • • = 
rifc = 2, the graph TB is isomorphic to the hypercube {0,1}^, in which three 
vertices are inserted in each edge, and the edge is duplicated into six directed edges 
as in Figure]^ The mapping between TB and the hypercube sends the tree 
to the point m, while the tree is interpreted as the vertex lying in the middle 
of the edge of the hypercube defined by the vector m' (this edge points in the Ath 
axial direction). 



Figure 4. Left: The hypercube {0,1}^; Right: the graph TB for 
k = 2 and rii = 71.2 = 2. 
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Let US now apply Theorem 3.5 to this example. For each I C [l..fc], let Wj 
be the strongly connected subset of B consisting of 0 and all vertices in the Tth 
petal of the bouquet for some i G I. It is easy to see that these sets are the 
only ones with nonzero multiplicity. By basic counting, it is immediate to see that 
m{Wi) = m{Wi,0) = ~ 1); since a spanning tree a of i? rooted at 0 is 

such that '0(a) = Wj if and only if the edge outgoing from the vertex i is (resp. is 
not) the one with smallest outgoing vertex for each i G I (resp. i ^ I). Moreover it 
follows from the interpretation in terms of rooted forests (Kirchoff’s theorem) that 
for i G [l..fc] one has 


( 6 . 1 ) 


detQwi = I Sj 

i zG[l..fe]\/ 


iG/ 3 = 1 


From Theorem 3.6 one thus obtains the value of the polynomial <i>B: 

n I E ^nlE^ ' 

/C[l..fc] \zG[l..fe]\/ / iG/ \i = l 

Equation( |2.4[ ) then implies that for any m G [l..ni] x [l..n 2 ] x ••• x [l..nfe] the 
generating polynomial of spanning trees oiTB rooted at a^ is given by: 


( 6 . 2 ) 


Z = 





Let us now examine more precisely the case ni = n 2 = ■■■ = rik = 2 and the 
link with the hypercube. Let Zm = Zm[y^,yl,ti,l < i < k) he the generating 
polynomial of spanning trees of the hypercube {0,1}^ rooted at m, where yl marks 
the number of edges in the tree mutating the i-th coordinate to the value j, and ti 
marks the number of edges of {0,1}^ that are parallel to the i-th axis and are not 
present in the tree. Then it is easy to see combinatorially (see Figure]^ again) that 
we have: 


(6.3) 


Z — ZfYii^SjXj , SiXj^ , “t“ Xj^ ). 


Therefore the value of the generating polynomial Tmiy^fiy}) = Z{y^^y]^ 1 ) can be 
recovered via the (invertible) change of variables x} + xf = 1, Sixj = 1 /°, and 


Sixf = y}, i.e. by substituting xj 


at 1 , and Si 
Vi+Vi 


y° +yl in ([O). 


Vi T 

v'i+Vi ’ 

We finally obtain the generating polynomial of spanning trees of the hypercube 
rooted at m: 


= n 


n i: 


.yO _J_ yl 

i=i /c[i..fc] \ie[i..k]\i 


(6,4) 


n»r‘ n 


\iej 

IJI>2 


in agreement with [31 Eq (13)] (see also [71 Thm. 3]). 

We note that a more refined enumeration can be obtained. First, let us now 
assign the weight w (instead of 1) to all the edges leaving the vertices vf, and let us 
replace the weights xf by wx^. Using Kirchoff’s theorem and a careful enumeration 








LAPLACIAN MATRICES AND SPANNING TREES OF TREE GRAPHS 


19 


of spanning forests of B, one can generalize (6.1) and prove that for I C the 
determinant det( 2 :/ — Qwi) is equal to: 




g wxi^ {w + z)'^'o ^ + z{w + z )"'0 

i=i 



(w + zY' 


This enables to apply Theorem 3.5 and obtain the full generating polynomial of 
forests of the graph TB. By extracting the top degree coefficient in w in the 
obtained formula, we obtain the generating function of spanning forests of TB in 
which roots can only be vertices “of type a”. In the case ni = n 2 = ■ ■ ■ = ni^ = 2, 
recalling that m(Wi) = 1 for all J, we obtain for this quantity the formula: 


I \ / iel j=l 

Now, the generating polynomial of directed forests on TB that have only roots of 
“type a”, and of spanning forests of the hypercube {0,1}^ are related combinato- 
rially by the same combinatorial change of variables as above, namely + ccf = 1, 
Six] = yY and Sixf = y}, that implies in particular that Si = y^ + y}. We thus 
obtain: 


Corollary 6.1 (0 Eq (3)]). The generating function of spanning oriented forests 
of the hypercube {0,1}^, with a weight z per root and a weight yf for each edge 
mutating the i-th coordinate to the value j is given by: 

Jc[l..fc] IGJ 


We conclude this section with a final comment. Of course, our proof of (6.4) or 


Corollary |6.1| via Theorem |3.6| is more complicated than a direct enumeration using 
Kirchoff’s theorem and an elementary identification of the eigenspaces. However, 
it sheds a new light on these formulas by placing them in the general context of 
tree graphs. Moreover, this places the problem of finding a combinatorial proof of 
these results and of our main theorem under the same roof. An indication of the 
difficulty of this problem is that as far as we know, and despite the progresses of [3] , 
no bijective proof of (6.4) (nor even (1.2|) is known. 
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